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We complete the formulation of a general framework for the analysis of high-order nonspherical 
perturbations of a four-dimensional spherical spacetime by including a gauge-invariant description 
of the perturbations. We present a general algorithm to construct these invariants and provide 
explicit formulas for the case of second-order metric perturbations. We show that the well-known 
problem of lack of invariance for the first-order perturbations with I = 0, 1 propagates to increasing 
values of I for perturbations of higher order, owing to mode coupling. We also discuss in which 
circumstances it is possible to construct the invariants. 
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I. INTRODUCTION 

The highly nonlinear character of the Einstein equa- 
tions of General Relativity requires the use of numerical 
and perturbative techniques to study problems of astro- 
physical interest. Processes like the merge of binary sys- 
tems (black holes and/or neutron stars) are promising 
candidates for gravitational wave detection in the next 
decade, and need careful and detailed simulation to ac- 
curately predict the waveforms that will be observed. 

The use of perturbation theory can be important both 
as a tool for simulating the evolution of a system in a 
state close to a known solution, or as a tool to extract 
the information of the amount and form of the gravita- 
tional radiation emitted. First-order perturbation theory 
is normally used in both cases, but it is possible to im- 
prove the accuracy of the results by using perturbation 
theory to higher orders. 

We have recently introduced a general framework to 
study generic high-order nonspherical perturbations of 
an arbitrary spherical spacetime [l[ . It can be considered 
as a generalization of the formalism of Gerlach and Sen- 
gupta (GS) for first-order perturbations around spherical 
symmetry |2| , in the sense that we use the same concepts 
and techniques: the background spacetime is decom- 
posed as the product of a 2-dimensional Lorentzian man- 
ifold and the 2-sphere; covariant notations are used on 
both submanifolds; and all perturbations are covariantly 
decomposed as series in tensor harmonics. However, the 
GS formalism also uses a gauge-invariant description, 
whereas our presentation in Ref. fij employed a fixed 
gauge (the Regge- Wheeler gauge |sj). The present arti- 
cle fills in this gap and completes the construction of the 
high-order generalization of the GS formalism by showing 
how one can construct gauge- invariant high-order pertur- 
bations. This allows one to treat the perturbations in 
any gauge while dealing with quantities and expressions 
that present the same form for all physically equivalent 
perturbations. 

The issue of perturbative gauge invariance is closely re- 
lated, though not equal, to the coordinate independence 
of General Relativity and has been dealt with in differ- 
ent ways already in first-order perturbation theory. The 



pioneering works by Sachs Moncrief Q, and Stew- 
art and Walker [6| established most of the basic ideas of 
the subject at first order, but offered different points of 
view about gauge invariance, what sometimes obscures 
the interpretation of the results and complicates any at- 
tempt of generalizing the analysis to higher orders. Brum 
and collaborators have shown [7[ that the geometrical ap- 
proach of Sachs [J], already rather constraining at first 
order, becomes even more restrictive at higher orders and 
hence only highly symmetric scenarios can be described 
in this approach. Here we will show that the alterna- 
tive approach by Moncrief [f| allows a greater and more 
convenient flexibility, so that we will adhere to it in or- 
der to construct metric gauge invariants at both first and 
second order around a spherical spacetime. 

This article, which can thus be considered a continua- 
tion of Ref. [1] (referred to as Paper I in the following), 
is organized as follows. Sec. |H] reviews the concept of 
gauge transformations and Sec. IIIII discusses the differ- 
ent notions of gauge invariance that have been employed 
in the literature, introducing a general method for the 
construction of gauge invariants. Sec. [IV] describes the 
notation adopted for spherically symmetric spacctimes 
and presents a brief summary on the structure of tensor 
harmonics (see Paper I for a more detailed description). 
Sec. [V] constructs gauge-invariant nonspherical pertur- 
bations around a spherical spacetime, using the method 
introduced in Sec. IIIII and is the central section of the pa- 
per. Finally, Sec. I VII contains our conclusions. Appendix 
lAlis devoted to the study of mode coupling and discusses 
in which cases the high-order gauge invariants can be con- 
sistently constructed. A particular case, already studied 
in the literature ;8| , is described in Appendix [Bj 



II. GAUGE FREEDOM IN PERTURBATION 
THEORY 

Perturbation theory in General Relativity considers a 
family of spacetimes (M £ , g s ) in which associated families 
of tensor fields fl £ are defined. The e = members of 
those families are referred to as the background members 
and will be denoted without a subscript, e.g. the back- 
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ground metric is g = go. All manifolds M e are assumed 
to be diffcomorphic. The main issue in perturbation the- 
ory is comparing tensor fields for a given nonzero value of 
e with their background counterparts. There exists dif- 
feomorphism invariance on each of the manifolds M s but, 
in addition, there is no preferred point-to-point identifi- 
cation mapping between any two such manifolds, so that 
the comparison of two tensor fields with different values 
of e is not an invariantly defined concept. This is the ori- 
in of the so-called gauge freedom in perturbation theory 



Let us call a gauge <f> e a family of point-to-point identi- 
fication diffeomorphisms from the background manifold 
to M £ : 



M 



Me- 



(1) 



Given a gauge cj) E we can now pull-back a generic tensor 
il £ on M s to a tensor <f)*£l e on M. This latter tensor 
can be compared with the background member f2 (at 
each point x € M), resulting in a (^-dependent concept 
of what a perturbation means. Assuming smooth depen- 
dence of all structures in e, we can define the perturbative 
expansion 
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where all terms of the equation are defined at the same 
point of the background manifold M. In particular the 
perturbations A?[f2] are tensor fields on M. The notation 
A<£ stresses the fact that, in general, it is not possible to 
define a perturbation without explicitly indicating which 
gauge 4> is used. For instance, the statement that a per- 
turbation vanishes is generically meaningless unless one 
specifies the gauge in which this occurs. Note also that 
the infinite series in Eq. $2$ arises from the simultaneous 
dependence on e of both fl £ and <p E . 

One then has to face the question of how the pertur- 
bations A^[0] vary under a change of gauge from cj> E to, 
let's say, ip e while keeping unaltered the family of tensors 
fl E . Such a gauge transformation will be described by a 
family Xe of diffeomorphisms on the background manifold 



M 



which clearly satisfy 



(3) 



(4) 



We emphasize that Xe is not a gauge, but a gauge trans- 
formation. 

The theory of expansions of families of general dif- 
feomorphisms on a given manifold has been developed 
in Ref. [3]. The most important result is that any 
one-parameter family of diffeomorphisms, in particular 
the gauge transformation % e , is equivalent to an infinite 
set of flows (a flow is a one-parameter group of diffco- 



morphisms) 
pressed in the following way 
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.}, so that Xe can be ex- 
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From the definition of the Lie derivative, it is straight- 
forward to see that the pull-back of any flow <f> £ acting 
on a generic tensor T can be expanded as 



F 



(6) 



where the vector £ is the generator of the flow. Making 
use of this formula and the decomposition ([5]), we can 
expand the right-hand side of Eq. ([¥]) in a power series, 
so that two gauge choices are related by 



00 °° °° £ k 1 +2k 2 +...+rnk m + ... 

«i=0fc2=0 k m =0 



fKl pk 2 r. 



(7) 

The vector fields { {2> £,..., fm} £, ...} are the corre- 
sponding generators of the flows { fl} $, f2} <£>,..., fm} <!>,...}. 
Substituting expansion {2J in the above equation, we get 
1 

n 

A$[0]-A3Pl]=X; 



' (n-m)\^ 2 k i...(m\) k ™k 1 \...k 
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where we have defined 



(k 1: ...,k m ) eN m ; Y, ih 



i=l 



Up to third order we obtain 
A^-A^fi] = £ {1}? r!, 



(9) 



A£[fi]-A£[n] = [c^ + ci^ 

+ 2£ {1}e A^[0], (10) 



(11) 



3(£ {2}e + £ 2 {1} J A [fi] 
3£ {1 > f A|[Q]. 



These formulas describe the effect of general gauge trans- 
formations on any high-order perturbation of a generic 
background tensor Q. They contain all the information 
needed to analyze the issue of gauge transformations in 
perturbation theory. 

In practical applications we can keep terms only up 
to a finite order e™, which projects the full group Q of 
gauge transformations into a truncated group { " } <? of 
rtth-order gauge transformations. Each of these trans- 
formations is described by a collection of n vector fields 
{ f f " } £}, an d we wm say that it is pure mth-order if 

all those vectors are zero except for {m} £. The important 
point is that, in general, composition of pure mth-order 
transformations is not pure mth-order, unless m = n. 
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For example, the composition of two generic second-order 
transformations described by { {1> £, (2} £} and { fl} £, {2} £} 
is described by the pair 

| {1 >£ m£ + (12) 

and hence the subset of pure first-order transformations 
{ {1} £,0} is not a subgroup of {2} Q. In fact, the group 
{1} Q is not equivalent to this subset, but only to a trun- 
cated form of it, and therefore it is important to dis- 
tinguish between first-order transformations { and 
pure first-order transformations { 0, ...0}. In general, 
the set 1 " } Q of all nth-order gauge transformations is a 
group, but the subset of all pure mth-order transforma- 
tions is not. The only exception is the reduced case 
rn = n of transformations of the form {0, ...,0, 
in which only a single linear term in Eq. ([5]) survives 
(this includes first-order perturbation theory as the case 
m = n = 1). There are more general subgroups of a 
given {n} G, like the subgroup of transformations of the 
form {0, ...,0, {m + 1} £, { ">£}, but they have less 
interest for our discussion and will not be considered in 
this work. 



III. THE NOTION OF GAUGE INVARIANCE 

Once we have defined the concept of gauge transforma- 
tion in Eq. ([3]), we discuss now the associated notion of 
gauge invariance of a family of tensors ft e under a group 
of gauge transformations. We will then find the inher- 
ited gauge invariance of the perturbations A^[ft] under 
the respective truncated version of that group. 

The most natural definition of gauge invariance was 
given by Sachs [4]: A tensor family ft e is identifica- 
tion gauge invariant (IGI) if the pull-back of its mem- 
bers to the background manifold is independent of the 
gauge, though the result still depends on e. That is, 
(f>*il E = ijj*£l E for all gauges (f> E ,ip E - This can also be 
interpreted as the invariance under the full group Q of 
gauge transformations. Perturbatively, a tensor family 
ft e is IGI up to order n if and only if A™ [ft] = A™ [ft] for 
all m < n and all gauges 4>, if) @- Again, this is equiva- 
lent to the requirement of invariance under the truncated 
group lni G defined above. This definition turns out to 
be too restrictive because, as it is well known 4, 6], only 
perturbations of vanishing tensors, constant scalars, or 
constant linear combinations of products of delta ten- 
sors can be IGI in first-order perturbation theory, since 
these are the only tensors with zero Lie derivative along 
every vector field. For higher, nth-order perturbations 
the problem becomes even worse because, apart from the 
background quantities, all of the mth-order perturbations 
with m < n must also be of the form that we have com- 
mented Q. In principle, this restricts to a very narrow 
physical scenario the possibilities that are left. 

Other forms of gauge invariance can be defined using 
subgroups of Q or For example invariance with 



respect to the reduced subgroup of pure nth-order trans- 
formations has been used in the past [l(|, by fixing the 
gauge perturbatively at all orders 1, n — 1, but not at 
order n. 

Only perturbations of highly symmetric backgrounds 
admit a complete description in terms of IGI vari- 
ables. Even at first order, significant limitations have 
been found: Stewart and Walker showed that, for vac- 
uum spacetimes, only backgrounds with Petrov type 
D are possible p, which fortunately includes the Kerr 
spacetime. In cosmology, only perturbations of static 
Friedman-Robertson- Walker (FRW) backgrounds can be 
described in terms of IGI variables pT| . For spherical 
backgrounds with matter, only first-order perturbations 
with axial polarity admit such a description [l2l . [l3[ , but 
not the complementary set of polar perturbations. This 
latter result is specially relevant for us, because we want 
to construct high-order gauge invariants of a spherical 
spacetime and in Paper I we saw that the polarities mix 
already at second order: hence there is no hope of getting 
a purely IGI description in general. Note that vacuum [5( 
and electro-vacuum [3] spacetimes with spherical sym- 
metry are very special cases (in particular included in 
the cited result for type D spacetimes), for which the 
programme of construction of gauge invariants can be 
further developed [la ]. 

On the other hand, when describing gravitational ra- 
diation in a vacuum, the Weyl tensor provides all the 
relevant geometrical information, and therefore many in- 
vestigations employ it as the basic object to be perturbed. 
Furthermore, the Weyl tensor defines a set of principal 
null directions, so that it becomes natural to decompose 
it using the Newman-Penrose formalism. The analy- 
sis of IGI is then simplified, but an additional type of 
gauge invariance is introduced, called tetrad gauge invari- 
ance, which requires invariance under (the 6-parametric 
Lorent z gr oup of) transformations among null tetrads 
0, EH, EH- We will not use this approach in this pa- 
per, but instead appeal to a different and more general 
notion of gauge invariance in which one makes use of 
an additional geometrical structure: a privileged gauge 
4> E . The basic idea is that, given a family of tensors ft e , 
one can select a privileged gauge to extract the physical 
information contained in this family and express this in- 
formation in terms of the pull-back of ft e in an arbitrary 
gauge. In other words, the gauge invariant is defined 
as the function(al) that provides the value of (f) E Q E in a 
generic gauge: 

£fi e = F [£fi e ] . (13) 

So, the gauge invariant is now supplied by the func- 
tional) F rather than by the family of tensors ft £ itself, 
as was the case for IGI. 

This notion of gauge invariants is similar to that of the 
constants of motion defined in Mechanics by the partic- 
ular values that the variables of the system take at some 
fixed instant of time [ItJ , or even to the notion of evolv- 
ing constants of motion recently introduced in Quantum 



4 



Gravity [l8| (although in that case one ought to con- 
sider a family of privileged gauges parameterized by a 
set of real numbers, rather than just one of them </> e ). 
In spite of appearing counterintuitive at first, this notion 
can be very useful in those cases in which the compu- 
tations can be carried out to completion, i.e. when one 
can obtain the explicit expression of the invariants in 
terms of gauge-dependent quantities, in our context, or 
in terms of time-dependent variables, in Mechanics. In 
other words, we need to determine the explicit form of 
the gauge transformation F — x* = (f)^ 1 o(f> £ for arbitrary 
<f> e . Whether this is possible or not essentially depends 
on the choice of gauge cf> £ . 

In practice, the privileged gauge is defined by imposing 
some conditions R £ on the pull-back 4> £ fl e of a particular 
tensor Cl £ . Therefore, 4> £ will be characterized as the 
gauge in which the tensor cf> £ tl £ satisfies some specific 
requirements. For this method to work satisfactorily, this 
privileged gauge choice has to be rigid. This means that 
the conditions R e [(f> e Cl s ] = must fix uniquely the gauge 
4> e , and so any further gauge transformation will violate 
those conditions. 

In perturbation theory, the invariants will then be the 
combinations '"^[{A™^]}] obtained by performing a 
gauge transformation from the perturbations defined on 
a generic gauge A™ [Q] to those defined in the rigid gauge 
A™[f2]. This kind of combination of perturbations have 
been characteristic of this approach to gauge invariance, 
starting with the pioneering work of Moncrief Q for 
nonspherical perturbations of Schwarzschild, where the 
Regge- Wheeler (RW) gauge was implicitly used as the 
privileged gauge. His work was later generalized by GS 
[2| to nonspherical perturbations of any spherical back- 
ground, also implicitly using the RW gauge. The same 
procedure has been employed by Bardeen Stewart 
11], and many other authors in their study of perturba- 
tions of FRW cosmologies. It can also be found in sev- 
eral recent investigations of second-order perturbations 
of vacuum d, or cosmological backgrounds [2(| • 

For instance, the first-order gauge invariants of a 
generic tensor f2 will be given by 



(14) 



where p^ is the vector field generating the first-order 
gauge transformation from <fi to (f>, so that this vector 
contains now information about our choice of privileged 
gauge 4>. Again, in practical applications the gauge </> 
is selected by imposing some rigid conditions R on the 
perturbations A^[fi] for some specific tensor f2, and such 
that no residual freedom is left in the choice of gauge: 



R 



= 0. 



In this way we get the equations 

rIa^Ci] + c p n] =o 



(15) 



(16) 



which must be solved for p M in terms of A^jTJ]. Sub- 
stituting the vector p^ obtained in this way, expressions 
(H3J) provide gauge invariants by construction. Note that 
when = f2, some of those expressions (or combinations 
of them) are trivial identities [equivalent to the require- 
ments (|16[) ]. This method for the determination of in- 
variants can be straightforwardly generalized to higher 
perturbative orders, as we will see in the following case. 

Since metric perturbations play a central role in our 
analysis, we choose the background metric as the 
tensor Q on which one imposes the conditions to fix the 
privileged gauge. We introduce the following compact 
notation for the perturbations of the metric: 



(17) 
(18) 



for a generic gauge <p and our privileged one </>, respec- 
tively. At first order we have that expressions for 
the metric become 



(19) 



The vector p^ is determined by demanding some condi- 
tions RlK.^] = which characterize the gauge <fi at first 
order. Then, the vector p^ is determined in terms of the 
components of by solving the equations 



R{h^ v + C p g^} = 0. 



(20) 



This completes the definition (fl9|) of the gauge invariant 
/C M „ as a function of h^,. 

Nonetheless we note that, owing to the presence of the 
Lie derivative, Eqs. (f2T)]) contain derivatives of the vector 
p M , so that their solution will involve in general integrals 
of the metric perturbations. Only when p^ can be de- 
termined explicitly in an amenable way from the metric 
perturbations we will have a useful form of gauge invari- 
ants. This fact will depend on the choice of the privileged 
gauge. In particular, we will see later that around spher- 
ical backgrounds the requirement of getting explicit and 
nonintegral expressions for the harmonic components of 
the vector p M will almost uniquely single out the RW 
gauge. We also point out that the same vector ob- 
tained by solving Eqs. ([2H)l . can now be employed to de- 
fine the gauge invariants associated with any other tensor 

as in Eq. (fT4]) . In addition, note that we can still in- 
terpret JC^ U as (the value of) the metric perturbations 
expressed in the rigid gauge <f> which satisfies conditions 

(HDD. 

At higher orders, and once a rigid gauge is chosen via 
some conditions {m} R for all m < n, one can obtain the 
nth-order metric invariants as 

f >/C„„ = { "V + C {nip g^ + in} n pv . (21) 

Since this equality reflects the effect of a gauge trans- 
formation, the source ^"^H^ is explicitly given by Eq. 
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© and depends on lower-order vectors {m >p^ and per- 
turbations {m ^h^ v with m < n, but not on f™>j/\ Be- 
sides, we remember that the source vanishes at first order 
( {iy Hfii/ = 0). On the other hand, the equation that one 
has to solve iteratively in order to determine the gauge 
vectors {ml p^, from m = 1 to m = n, takes now the 
expression 



{mi R[ {mi h^+C im}p g^+ {m} ?V] =0. 



(22) 



In particular, when all the conditions {m) i? have the same 
linear functional dependence on their arguments (for in- 
stance because they arise from the perturbative expan- 
sion of just one set of exact linear gauge conditions on 
the metric), Eq. (|2"2"|) will have the form (|2T)| but with 
the source term {m} R[ {m} Ti.fj, u ]. Therefore, the solutions 
of these equations will be constructed essentially in the 
same way. 

Nakamura has suggested a similar approach [2l| to con- 
struct high-order gauge invariants. He starts from the 
basic assumption that a splitting equivalent to Eq. (|21|) 
is given from the outset, separating the metric pertur- 
bation {n) h^ v into its gauge-invariant part {,l} /C M „ and 
gauge- variant part (containing the vectors {m} p^), with 
the vectors {ml p^ satisfying some set of requirements. No 
proposal is made, however, on how such a splitting can 
be attained. Our scheme goes beyond that proposal, giv- 
ing a constructive and general prescription to generate 
the vectors {m) p^ from the choice of a rigid gauge, in 
such a way that the requirements imposed on t" 1 ^ are 
automatically fulfilled. 

After determining the vectors { lli p^(h), {n) p^(/i)}, 
the perturbations of any tensor field, and in particular 
those of the stress-energy tensor {n} ^ can be taken to 
its gauge-invariant form '"'T^^ just by applying a gauge 
transformation parameterized by the above vectors: 



n 



E 



i 



E ( n _ m )| 2 fc 2...(m!) fc ™/c 1 !...fc m ! 

rn—0 x ' (Km) 



t-'inp'"*-' {m} p M^' 



(23) 



In this way we will get a tensor {n} T /J!/ (^, h) whose 
dependence on the perturbations {ml, I , M „ and {m} h^ v 
(m < n) will not change when any gauge transformation 
is applied to them. In Sec. |V]we will use these techniques 
to compute the metric and matter gauge invariants for 
perturbations of a spherical background spacctime. 



+ r z (x u ) lah {x d )dx a dx\ (24) 
T^{x D .x^dx^dx" = t AB (x D )dx A dx B + ^r 2 (x D ) 

x Q{x D ) lab {x d )dx a dx\ (25) 

where Greek indices, capital Latin indices, and lowercase 
Latin indices take values in the ranges {0, 1, 2, 3}, {0, 1}, 
and {2,3}, respectively. In this way, qab is the metric 
of the manifold A4 2 , whereas j a b is the round metric on 
the sphere, and r is a scalar function defined on A 7 ! 2 . In 
order to avoid working explicitly with logarithms of the 
function r, it is usual to define 



va 



r.A 



= (hir) jj4 . 



(26) 



We will use the following notation for the covariant 
derivatives associated with the different introduced met- 



rics: 



9fiu;p = 0, 9AB\D = 0, 7ah:d = 0. (27) 

In Paper I we generalized the so-called Regge-Wheeler- 
Zerilli basis of harmonics [3, [12] to any rank. Specifically, 
we introduced the following two independent rank s sym- 
metric and traceless tensors on the sphere 



J Z a±...a s 



X, 



I oi 



(^ m :ni...aJ 

b rym 

e (ai ba 



STF 



.a s ) ! 



(28) 
(29) 



where Y t m are the scalar harmonics and the superscript 
STF means the symmetric and tracefree part. This defi- 
nition is valid for —l<m<l and 1 < s < I. In all other 
cases, these harmonics are defined to be identically zero, 
except for s = 0, when Z™ = Y, m . These tensor harmon- 
ics provide two families with distinct polarity: ZJ n ai ,.. aB 
is polar, whereas A ; m ai ... Qs is axial. Under a parity trans- 
formation, the polar and axial harmonics change sign as 
(—1)' and (— respectively. 

A basis for tensors of rank s on the sphere is then 
formed by the two tensors {ZJ n ai ,,, as , XJ n ai ,,, as } and 
some linear combinations of products between a basis of 
tensors of rank (s — 2), on the one hand, and the metric 
jab and the Levi-Civita tensor on the sphere e a b, on the 
other hand. 



IV. SPHERICALLY SYMMETRIC SPACETIMES 



GAUGE INVARIANT PERTURBATIONS OF 
A SPHERICAL BACKGROUND 



In the following we will restrict our discussion to spher- 
ically symmetric backgrounds and employ the notation 
introduced by GS @|. We decompose the background 
manifold M. as the product M 2 x S 2 , where A^ 2 is a 
two-dimensional Lorentzian manifold and S 2 is the two- 
sphere. Without loss of generality, any spherically sym- 
metric metric and stress-energy tensor can be written in 
the block-diagonal form 0: 

g tliy (x D ,x d )dx'"dx u = g A B{x D )dx A dx B 



A. Harmonic decompositions 

When working on spherical backgrounds, it is conve- 
nient to decompose our geometrical objects in series of 
spherical harmonics in order to remove all angle depen- 
dencies from the equations of the problem. We therefore 
adopt from now on the following decomposition for the 
different tensors of interest: 



6 



,tj UV 



1 tuy 



E 

l, in 

E 

l,m 

E 

l, in 

E 

l, in 



H l AB 

Sym. 
Sym. 
Sym. 



(30) 
(31) 

(32) 
(33) 



r 



The notation "Sym." indicates that the considered ma- 
trices are symmetric. In the most general case, the label 
I runs over all nonnegative integers, whereas m has the 
usual restriction —l<m<l. Note also that we have 
adopted the convention that polar (axial) harmonic co- 
efficients are denoted by capital (lowercase) letters. 



B. Privileged gauge 

Using the previous harmonic expansions we can de- 
compose the nth-order gauge transformation (|2ip from 
the perturbations of the metric in a generic gauge {n ^h^ v 
to the perturbations { " } /C M1/ in some privileged gauge, 
that we still have to specify: 



I > : 



} HPAB + {ni Pl n A\B + in} Pl n B\ 



+ {ny H?AB 



' <H AB, 

{ n }jC m _ {ny K m + 2v Alni P l m A -l{l + l) lni P{ 



(34) 
(35) 



I > 1 



/ > 2 

{n 



Ta = M Hr A + M pr A + r 2 M P[ n \A 

+ {n} n?Ai 

k?a = {n} h™ a + r 2 { " W + {n} h™A, 



•m 



(36) 
(37) 

(38) 
(39) 



The cases I = 0, 1 are special inasmuch as some of the 
tensor spherical harmonics vanish for these values of I, 
and hence some of the relations that one obtains for I > 2 
are trivially void in these cases. 

In general, the privileged gauge at order n will be de- 
termined by four conditions {n} R, imposed on the ten 
components of { ™ } /C M „. These conditions must rigidly 
fix the gauge vector (j33|) as a function(al) of the generic 
metric perturbation {n) h^ v . In order to obtain amenable 
gauge invariants, we would like to select a gauge such 
that this functional dependence is local with respect to 



the Lorentzian manifold A4 2 (whereas the dependence on 
S 2 is determined by means of the expansion in tensor har- 
monics). After replacing the obtained gauge vector, the 
other six components of {n} /C^„ will provide then local (in 
M 2 ) gauge invariants by construction. For simplicity, we 
will employ from now on the adjective "local" exclusively 
in the sense of locality in A4 2 , so that no integration over 
this manifold is involved. 

From the above expressions, it is not difficult to see 
that there are essentially two procedures to end with a 
vector whose components present a polynomic de- 

pendence on those of {n} K,^ u and '""'/i^^, for m < n. In 
both cases, the components f™>p™ and are found 

solving Eqs. (j3"8)) and (f3"9"| . respectively. But for the vec- 
tor component {n} P[ n A one has two options. One can 
fully determine it from Eq. (|36p . or one can instead ob- 
tain the projection v A {n} P™ a from Eq. (|35p and the 
remaining projection t A f™>P™ A from Eq. (|3"6"|) , with t A 
being any vector transverse to v A . 

This fact suggests two possible sets of conditions to 
select our privileged gauge, both of them imposed on the 
tensor {n} /C^ v . The first one consists of the following 
requirements 



1>1: { ' 
l>2: fr 



'/C 



l A 



0. 



} icp = o, 



In} m 

K l 



0. 



(40) 
(41) 



This is the RW gauge used at first order @ . This gauge 
has been extensively employed in the study of perturba- 
tions of spherical backgrounds and leads to a full metric 
whose pullback satisfies |l| 



{(j> e gs)Ab:c.g bc = 0, ((f) e g £ )ab = K e g a b, 



(42) 



for some generic scalar field K e defined in the background 
four-dimensional manifold Ai . It is well known that this 
gauge is not completely rigid because it does not impose 
any restriction on some of the I = 0, 1 components of the 
metric perturbations, at any order. For the rest of this 
subsection, we will restrict our analysis of gauge invari- 
ance to transformations whose generators do not contain 
the corresponding 1 — 0,1 modes (both at order n and 
lower orders). We defer the analysis of the general case 
to the next subsection. 
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With the choice of the RW gauge, Eqs. (J36J) , (38]), and 
([55)) are the conditions that must be solved to obtain the 
gauge vector { ">p' 1 , which turns out to be 



I A 



12 rt i 



0}pm _ _ 



I A, 



/ {"}/,"! 



i^ H rn A 

(43) 

(44) 
(45) 



Substituting these expressions in Eqs. (|34|. (|35j) . and 
Q37[) we obtain explicitly the high-order metric gauge in- 
variants. 

The axial part of the second algebraic gauge is equiva- 
lent to the RW gauge, and therefore { ™>g ; m is given again 
by Eq. (|4"5"|) . As a consequence, the axial invariant is the 
same as before, namely { " } k™,4- However, the polar part 
is now determined by the following conditions: 



I > 
I > 1 
I > 2 



M K7 = 0, 



{n 



o, 



(46) 
(47) 
(48) 



which generalize the Campolattaro-Thorne gauge K = 
[23l | for the I — 1 component of first-order perturbations. 
The exact gauge would be 



{4> e 9e)Ab:c9 bC = e A Bt B H £ 



{<t>s9s)ab = 9ab, (49) 



for some scalar field H £ in the background M. . The above 
comments about the lack of rigidity are also applicable 
to this gauge. That is, the gauge generators are assumed 
to have no component with harmonic label I = 0, 1. 

For the polar components of the vector {ret p^, we solve 
conditions (|46|) and obtain 



} A P m „ 

J "l A 



t A P? 



PT 



} G7 



1 



-t" 



A in 



1 



{n 



2 



\A 



(50) 

(51) 
(52) 



Substituting the above relations in the respective defi- 
nitions (|34|) and (|36|) . one arrives at the corresponding 
polar invariants K,™ ab and v A IC] n A- 

In conclusion, we have seen that there are two choices 
of gauge that lead to a polynomic and local expression 
of the gauge vector {ni p^ in terms of the metric pertur- 
bation {7l} /i^„ in an arbitrary gauge. For both of them, 
one obtains a set of local gauge invariants, under the re- 
striction I > 2 on the modes of the gauge generators, at 



all perturbative orders. These two sets are inequivalent. 
Nevertheless, the RW gauge does not need the introduc- 
tion by hand of a vector field t a transverse to va , so that 
one can consider it to be more natural, at least as far 
as this kind of geometric considerations are concerned. 
Because of these reasons, we choose the RW gauge as 
our privileged gauge from now on. In Paper I we already 
computed the evolution equations for the perturbations 
in the RW gauge. Of course, those equations are also 
valid for the gauge invariants that are tied to this gauge. 

The fact that there exist privileged gauges that pro- 
vide local gauge invariants is related to the spherical 
symmetry of the background, which allows the choice 
of an adapted system of coordinates in which the metric 
has a block-diagonal form. Besides, Kodama, Ishibashi, 
and Seto [13] have shown that it is possible to generalize 
the first-order GS invariants for any metric which can be 
written in block-diagonal form in an (m + n)-dimensional 
manifold M. m x K n with a maximally symmetric block 
in K, n . This means that a higher-dimensional form of the 
RW gauge is always attainable, though the corresponding 
conditions on the metric perturbations may be different. 



C. The particular cases / = 0, 1 

It is well-known that, at first order, one cannot con- 
struct local gauge invariants for / = 0, 1 by the methods 
explained in the previous subsection. This is because 
some of the equations (|34| - [39|) are not present in those 
cases, and therefore it is not possible to attain a local 
expression for the components of the gauge vector 
in terms of those of the metric perturbations. Of course, 
gauge conditions different from the RW ones can be im- 
posed on the metric perturbations and hence one can 
obtain from them the associated gauge invariants, but 
these will be nonlocal because the gauge vector will be 
given by an integral expression (over A4 2 ) of the metric 
perturbations. Whether this is useful or not will depend 
on the particular application that one is studying. 

The same obstruction appears as well at every order in 
perturbation theory, so that one cannot get local gauge 
invariants for I = 0, 1 at any order. However, mode cou- 
pling makes the problem worse: the existence of lower- 
order modes with / = 0, 1 may prevent the construction 
of higher-order local gauge invariants with I > 2. This 
will happen when such lower-order modes have a nonzero 
contribution to the sources TC of Eqs. (|3~4"H3"9")) . 

The restriction on the gauge generators introduced in 
the previous subsection is a valid alternative for first- 
order perturbation theory: the GS gauge invariants re- 
main unchanged under the restricted group fl} <?° , where 
the diamond indicates that no first-order generator with 
I = 0, 1 is included. However, removing all I = 0, 1 gen- 
erators is already inconsistent at second order if we de- 
mand invariance under a group of transformations, be- 
cause the I = 0, 1 components of the vectors {2> £ will be 
unavoidably generated by coupling of first-order gauge 
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modes [cf. composition (p~2|) ] . Fortunately, those offend- 
ing gauge modes act only on the Z = 0, 1 second-order 
perturbations, for which invariants cannot be constructed 
anyway. All other second-order perturbations admit a 
gauge-invariant form, as given in the previous subsec- 
tion, under the gauge group {2} CP where again the dia- 
mond denotes that no first-order I = 0, 1 gauge mode is 
included, but all other first and second-order (including 
I = 0, 1) modes are allowed. That is one of the main re- 
sults of this investigation: gauge-invariant perturbations 
can be consistently and simultaneously constructed for 
all Z > 2 modes at second order, all being invariant under 
the group {2} Q°. 

The situation at third and higher orders is more re- 
strictive. As a summary, working at order n > 2 only 
a finite number of gauge generators can be included in 
the invariance group at orders 1, ...,n — 2. The unavoid- 
able presence of second-order gauge modes with I = 0, 1 
couples to any metric perturbation at order m and with 
harmonic label Z, preventing the construction of a gauge 
invariant of order m + 2 and label Z or Z ± 1, with respect 
to those gauge modes. The detailed analysis is shown in 
Appendix 



D. First-order gauge invariants 

We now particularize our discussion to first-order per- 
turbations and reproduce the GS first-order gauge invari- 
ants , which are a generalization of the invariants in- 
troduced by Moncrief [5|] for a Schwarzschild background. 

Since the source {1} H M „ vanishes, from Eqs. (f3"4"| . (pj)) . 
and (|37p one gets the metric invariants 

{1) ICr A B = {1) HP AB + ^Pr A \B + {11 P/™B|A,(53) 

- 1(1 + 1) {1} P™, (54) 
^VA + r^^rV, (55) 



and (|37|) for n — 2, whereas the matter invariants, con- 
structed from the perturbations '"^^ = A"^], are 
encoded in the tensor {2} T M ^ defined by 

+ 2£ {1}p fl >^. (59) 

For the sake of brevity we will only provide explicitly 
the form of the metric gauge invariants, but not that of 
the matter invariants, which can be computed along the 
same lines. 

Our aim is to compute the harmonic coefficients 
of the source 
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= C\ 1) 9 ^ + 2C {1}p {1 V- (60) 



These coefficients will have the general form (we obviate 
the n — 2 label for simplicity): 



= E E ^ 



1,1 m,m 



E E 



mmm 

i i i A ' 



E E 



mmm 
Til' 



n_/m 

rti AB 

n_jm 
H-l A 

HT 

Ka = -EE <- ^ 

hT 



l l l AB, 



1,1 



E E ^ 



mmm 
Til' 



mmm 
l l l A, 



U m,m 



*EE^ 



mmm 



(61) 
(62) 
(63) 
(64) 
(65) 
(66) 



K l A 



1,1 "i,m 



where 



{Dp 



I A 



{1} G 



l \A, 



1l 



2 
1 



2 n l ■ 



(56) 

(57) 
(58) 



Using these expressions for the gauge vector, it is possible 
to construct the first-order gauge invariants correspond- 
ing to any other (matter or metric) tensorial object. 



E. Second-order gauge invariants 

The second-order metric invariants can be obtained 
by substituting expressions (|!3Tt4"51) in Eqs. (JMJ), (j3"5)) . 



with Z and Z being the harmonic labels corresponding to 
the two first-order modes that get coupled. They are 
independent and run over all nonnegative integers with 
the restriction |Z — Z| < Z < |Z + l\- In addition, one 
must impose the usual restrictions on the values of the 
labels rh and rh. We have also defined the alternating 
sign e = (— 1)' + '~ ; . As it is shown in Appendix IA1 this 
alternating sign encodes the mixing of polarities caused 
by the product of tensor harmonics in higher-order per- 
turbation theory. For a detailed explanation of this effect 
and the use of the sign e we refer the reader to Paper I. 
As in that reference, all computations have been carried 
out with the computer algebra system xTensor [2q ]. 

From now on and without loss of generality, we will 
suppose that there are only two first-order modes that 
will be denoted as (IC,p) and (IC,p). The tensorial sources 
are given by 
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l l I l AB 



J I AB 



-AE-% [Pt AB + 2q l(B k A) - {PP (B )\ A) 



r z P\ A P\B 



2E Opil \ 2pC \(A^B)C 



P C K 



AB\C 



P C \(bPa)\c 



P C \aPc\b 



r q\AQ\B 

P C P, 



(A\B)C 



iqK AB + 2P\ [A k B) - 2r 2 q\ (A P\ B) + qP {A \ B) + q\{ A P B )\ 
I 



(67) 
(68) 



The real coefficients E come from the product of har- we have defined the normalization factor 
monies 1] that is implicit in Eq. (|60p. They are defined 
as 



lm _fc(ZM*W,l«l) 



tjs l m 



k(l",\s + s'\) 



\ \) f-im' mm' +m r~is' s s' -\-s i aci\ 

Th - °/' i i" °r 1 1" ' \ 0J ) 



where c^™ 2 ™ 1+m2 are Clebsch-Gordan coefficients and 



k(l,s) 



'(2l + l)(l + s)\ 



2 s + 2 tt(/-s)!' 
The vectorial sources are decomposed as 



(70) 



n i i i a 



rp2 lm 
^-liml 



-2(qk A + qk A ) + (P A P + P A P) + r 2 {3PP ]A + PP\ A ) + 2q(r 2 q)\ A + r 



6 (QQ 



-E° 1 \^2l(l + l)(qk A - qk A ) + AP B iC AB + 4r 2 P| A £ 

1(1 + 1) (PP A + PP A + 3r 2 PP\ A + r 2 PP {A + r 2 qq ]A - r 2 qq ]A 



- 2r 2 



(p B Pib) 



\A 



Av B P B P\ A 



- 2 (P B P A \ B + 2P B P BlA + P B P MB ) j, 



i i i A 



-iE 2 _ x % {2 (Pk A - Pkj^j + qP A ~ qP A + r 2 qP\ A + 3r 2 Pq ]A - 3r 2 qP [A - r 2 Pq\ A } 
iE\% [2r 2 tq ]A - 1(1 + 1) (Pk A + Pk A ) + 2 (k B P B \ A + P B R A y 



+ ^-J^ (r 2 qP\ A - r 2 qP\ A + 3r 2 Pq\ A + r 2 Pq\ A - qP A + qP A 



r 2 (P B q\ BA + P B \ A q\ B + ±v B ' P B q\ A ) } 



and finally the four scalar sources are given by 

(+)H rnmm = _ AE ^% (qq + Pp) + E^ 1 \% t - -I P/v' 



J 2 frhl 



jrhl 



(71) 



(72) 



1(1 + 1) + 1(1 + 1) PP + ^P A (r 2 P) lA + \P A P A 



+ E^{ - 2 [([ + l)re + ^(^) |A + ^^(r 4 F) | , 



- 1(1 + 1)1(1 + l)PP - 2P A [(P B v B )\ A + 2P B v A v B ] j, 
(-mffT = %%E^ 2 %Pq-2iE^% [\p A k A - 2qt - 1(1 + l)Pq + \p A (r 2 q)\ A 



wfL mmm = 2r 2 E^% (qq + Pp) + v 2 E\ % { (/ + 2) (1 1) 



PP - qq 



- 2P A P {A - —P A P A 



{•■ 



2r 2 Ef i t l \2PiC + 1(1 + 1) 



qq + 2PP 



( — )n_/rnrn')n 

1 1 1 



2ir 2 E, 



J 3 frhl 

iE\%{AP A k A + (i+2){i-l)r 2 



2Pq - Pq 



P A (r i q) lA 



Pq + Pq 



2r 2 P A q\ A } . 



(73) 
(74) 

(75) 
(76) 
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These formulas allow us to construct the gauge-invariant 
form (under the restricted group {2} (/° of transforma- 
tions) of any second-order metric perturbation in a co- 
variant way. Note that we have not restricted the har- 
monic labels corresponding to the first-order metric per- 
turbations. In particular we allow the values I = and 
I = 1 for them, although then some of the harmonic co- 
efficients in the formulas will vanish. 

Appendix [B] studies the particularization of these 
gauge invariants to the case of second-order perturba- 
tions of a vacuum in Schwarzschild coordinates when one 
has at first order a single polar mode with 1 = 2, com- 
paring the results with those of Ref. [|| . 



VI. CONCLUSIONS 

We have completed the construction of a framework 
for the analysis of high-order gauge-invariant nonspheri- 
cal perturbations of a general spherical spacetime, with 
any type of matter content. Our work started in Pa- 
per I by generalizing to higher orders the GS formal- 
ism for first-order perturbations, although the perturba- 
tions were described in that work in a particular gauge. 
In this article we have given explicit algorithms to con- 
struct gauge-invariant perturbations at any order, focus- 
ing again on the second order, for which we have provided 
full expressions of the metric gauge invariants. 

Before doing that, it has been necessary to discuss in 
detail what we mean by gauge invariance. The simplest 
or most natural form of gauge invariance (introduced by 
Sachs [H) turns out to be too restrictive in spherical 
symmetry, even more when dealing with high-order per- 
turbations. We have used an alternative form of gauge 
invariance, based on the work by Moncrief [H on per- 
turbations about Schwarzschild. Borrowing ideas from 
Bardeen [li| , we interpret Moncrief 's gauge invariants as 
being associated with a choice of preferred gauge. This 
reinterpretation has allowed us to extend the construc- 
tion of gauge invariants to higher orders in perturbation 
theory, our main goal in this work. 

The question then arises of what has been really gained 
in defining gauge invariants if these are associated with a 
preferred gauge. First of all, it is important to emphasize 
that these are true invariants, in the sense that they do 
not change under gauge transformations of the perturba- 
tions, and therefore in a given problem we do not have to 
worry about the gauge in which the perturbations actu- 
ally are described. Secondly, the freedom to choose the 
preferred gauge is severely restricted in practice by the re- 
quirement of obtaining amenable explicit expressions for 
the invariants. For example in spherical symmetry the 
demand of arriving at explicit expressions that are local 
with respect to the dependence on JM 2 almost singles 
out the RW gauge, which essentially renders the associ- 
ated invariants unique. For nonsymmetric background 
spacetimes this type of gauge invariants, with an explicit 
and manageable form in terms of the perturbations in a 



generic gauge, will generically not exist. 

Unfortunately the RW gauge cannot be imposed on 
perturbations with I = or I = 1 at all orders, and this 
obstacle propagates to higher values of I through mode 
coupling. In this way, the invariants associated with the 
RW gauge remain unaltered under the set of gauge trans- 
formations whose generators exclude all the modes with 
harmonic label I equal to or 1. In general, however, 
this set of transformations does not have group structure 
owing to mode coupling. Nevertheless, at second order, 
the group {2i Q° (which does not contain the I = 0, 1 
first-order gauge generators) allows the construction of 
all gauge invariants with I > 2 both at first and second 
orders, invariants that are exclusively associated with the 
RW gauge and whose expression we have been able to de- 
duce for the most general case (i.e., for any matter model 
and any background spacetime). 

Our framework is now ready to be applied to specific 
matter models and we are currently working on its ap- 
plication to the study of second-order perturbations of a 
time-dependent ball of perfect fluid, with a general (2- 
parameter) equation of state. 
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APPENDIX A: MODE COUPLING AMONG 
GAUGE GENERATORS 

We have constructed second-order gauge invariants un- 
der the group {2 >C?° starting just from the non(fully-)rigid 
RW gauge. This is not possible at higher orders and actu- 
ally one has to restrict to a finite set of lower-order modes 
both in the gauge generators and in the perturbations in 
order to define some form of gauge invariance. This is 
because the presence at first order of any gauge mode 
I > 2 will generate, just by self-coupling, the second- 
order modes with harmonic labels and 1. It will then 
be imposible to construct the gauge-invariant form of a 
third-order perturbation whose source TL contains a term 
coupling any of these second-order I = 0, 1 modes with 
any first-order mode. But those sources generically con- 
tain all possible couplings, and so only a problem in which 
we restrict the number of first-order gauge modes allows 
some form of gauge-invariance at third order. In this 
appendix we will analyze generic mode coupling around 
spherical symmetry, starting at second order and then 
proceeding to higher orders. We will later give some 
bounds on the number of modes that can be present at 
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first order to allow for the construction of a nth-order 
mode with label /. 

The second-order Z-mode will get a contribution from 
a pair of first-order modes I and I if two conditions are 
obeyed. On the one hand, the harmonic labels must be 
related by the usual composition formula 

\l-l\<l <l + l. (Al) 

On the other hand, mode coupling must conserve parity. 
To any harmonic coefficient with label I, we associate 
a polarity sign a such that, under parity, the harmonic 
changes by a sign u{— 1) . Polar (axial) harmonics have 
cr = +l (cr = — 1). Then, parity conservation implies the 
second condition: 

{-l) T+t - 1 = e = aaa, (A2) 

where we have made use of the alternating sign e. There 
is a special case in which the coupling of two modes sat- 
isfying Eqs. (|A1|) and (|A2|) does not contribute to a 
second-order mode, and the reason comes from the prop- 
erties of the Clebsch-Gordan coefficients that appear in 
the product formula for the tensor harmonics (see Paper 
I). In axisymmetry (fh = m = 0) the Clebsch-Gordan 
coefficients, and as a consequence the E-coefficients that 
couple the modes [defined in Eq. (fBU)) ]. vanish if I + I + I 
is odd. A highly geometric derivation of the perturbative 
structure under the assumption of axisymmetry is given 
inRef. 

This analysis can be extended to higher orders. In 
particular, the parity condition will be that a collection 
of k modes with harmonic labels {Zi, ...,/*} and polari- 
ties {(Ti, Ofc} will contribute to the mode (I, a) only if 

(-i)V = nt 1 (-i)^ i . 

Let us finally consider the case in which we have a 
first-order finite collection of modes with their harmonic 
labels taking all the values from I = 2 to I — Z max , with 
contributions from both the polar and the axial sectors. 
Coupling of these modes at order n will generate some 
new modes, following the above rules, so that the high- 
est value of their harmonic label will be n/ max . The con- 
struction of the gauge invariants under the correspond- 
ing group of transformations and tied to the RW gauge 
is only guaranteed for those modes with harmonic label 
greater than (n — 2)Z max + 1. This number comes from 
the coupling of the (n — 2)th-order (n — 2)Z max -mode with 
the second-order I = 1 mode. 



APPENDIX B: COMPARISON WITH THE 
RESULTS OF GARAT AND PRICE 

Carat and Price calculated the contribution of two ax- 
isymmetric (to = rh = 0) first-order polar quadrupolcs 
(I = I = 2) to the I = 2 multipole of the second-order 
Zerilli gauge invariant [8[. Their calculation was done 
in Schwarzschild (t, r) coordinates. In order to translate 
their coordinate-dependent expressions to our covariant 
notation, we introduce the following frame of vectors: 

*-r%)\ (bd 



where / = yr\ A r\ A = \/l — 2M/r and M is the back- 
ground black-hole mass. They are normalized in the fol- 
lowing way, 

r A r A = 1, t A t A = -1, r A t A = 0. (B2) 

The metric and Levi-Civita tensor are given by 

9ab = ~t A t B + r A r B , (B3) 
£AB = r A t B - r B t A . (B4) 

Note that the relation between our frame vectors and the 
coordinate 1-forms is 

ta = r\ A , t A = -ft ]A . (B5) 

Another question to consider in order to compare our 
results with those of Carat and Price is the definition of 
the tensor spherical harmonics. Firstly, their normaliza- 
tion is different, which makes their harmonic coefficients 
to be a factor a -1 = 2-^/^/5 greater than ours. Secondly, 
there is a mix of harmonic coefficients G and K because 
Garat and Price employed a harmonic proportional to 
Y§ . ab instead of using Z^ab- 

The translation from the harmonic coefficients used 
by Garat and Price (which appear below in the left-hand 
side and in their notation) to our coefficients is 



o-(n) 
H ' 


-» at A t B{n} H AB , 


(B6) 


R (n) 


- ar A t B ^H AB , 


(B7) 


ffW - 


- ar A r B ^H AB , 


(B8) 




-» aft A ^H A , 


(B9) 


/4 n) - 


_> af-^^HA, 


(BIO) 


K {n) . 


-» a(3 f " } G+ ,n} K), 


(Bll) 


G (n) _ 


-» a M G. 


(B12) 



In our notation, the n-th order Zerilli variable, in terms 
of the GS invariants, is given by 

<">IT Z = -J^{fr B ^K, AB -r^1C ]A )r A 

+ r w /C. (B13) 

In particular, at second order we could add to this 
function any quadratic combination of first-order modes, 
which would simply change the source term of the Zerilli 
equation. Nevertheless, this is the simplest Zerilli vari- 
able one can define. Using the sources and ([7T1 - 
[76)) we have computed the second-order GS invariants, 
compared our Zerilli function ()B13|) with that found by 
Garat and Price, and verified their formula (41) of Ref. 
1- 
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